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ABSTRACT 

We present a simple model that accurately describes various statistical properties of 
peculiar velocities of dark matter halos. We pay particular attention to the following 
two effects. First, the evolution of the halo peculiar velocity depends on the local 
matter density, instead of the global density. Secondly, dark matter halos are biased 
tracers of the underlying mass distribution, thus halos tend to be located preferentially 
at high density regions. For the former, we develop an empirical model calibrated with 
A^-body simulations, while for the latter, we use a conventional halo bias model based 
on the extended Press-Schechter model combined with an empirical log-normal prob- 
ability distribution function of the mass density distribution. We find that, compared 
with linear theory, the present model significantly improves the accuracy of predic- 
tions of statistical properties of the halo peculiar velocity field including the velocity 
dispersion, the probability distribution function, and the pairwise velocity dispersion 
at large separations. Thus, our model predictions may be useful in analyzing future 
observations of the peculiar velocities of galaxy clusters. 

Key words: galaxies: clusters: general - cosmology: theory - dark matter - large- 
scale structure of universe 



' 1 INTRODUCTION 

X 



In the standard structure formation scenario, dark matter 
halos are thought to be formed via the gravitational amplifi- 
cation of initial small density fluctuations, and their motion 
results from gravitational forces acting over large scales. As 
a consequence, the root-mean-square (rms) peculiar veloci- 
ties of halos at the present time depend on the total amount 
of matter in the universe, S7o, and on the normalization of 
the matter power spectrum, as, as oc cr^o' 6 - In addition, 
linear theory predicts that the halo peculiar velocity evolves 
as oc a(t)dD(t)/dt where a(t) is the scale factor normalized 
unity at the present time and D(t) is the linear growth rate 
(see Section 2) . Thus, a measurement of the rms peculiar ve- 
locities of halos can be used as a test of structure formation 
scenarios in principle. 

A measurement of the peculiar velocity of a dark mat- 



* This paper was published in Mon. Not. R. Astron. Soc. 343, 
1312-1318 (2003). Owing to an error in numerical computations, 
some incorrect results were presented there. Erratum is to be 
published in Mon. Not. R. Astron. Soc. Conclusions of the original 
version are unaffected by the correction. This version supersedes 
the original version. 



ter halo is possible through its baryon contents. Because a 
cluster of galaxies is formed in the gravitational potential of 
a massive dark matter halo, its peculiar motion traces that 
of the hosting halo. Therefore, the rms peculiar velocity of 
clusters of galaxies provides a reasonable measure of that of 
halos. So far, empirical distance indicators such as the Tully- 
Fisher and the D n -a relations have been used to measure 
the peculiar velocity of relatively nearby systems (e.g. Bah- 
call, Gramann & Cen 1994; Lauer & Postman 1994; Bahcall 
& Oh 1996; Moscardini et al. 1996; Borgani et al. 1997; 2000; 
Watkins 1997; Dale et al. 1999; Hudson et al. 1999; Colless 
et al. 2001). Future observations of the kinematic Sunyaev- 
Zel'dovich effect will provide measurements of line-of-sight 
peculiar velocity for distant systems (Sunyaev & Zel'dovich 
1980; Rephaeli & Lahav 1991; Haehnelt & Tegmark 1996; 
Kashlinsky & Atrio-Barandela 2000; Aghanim, Gorski & 
Puget 2001). 

Theoretical predictions of the rms peculiar velocity of 
halos have been made by applying the peak theory (Bardeen 
et al. 1986), in which halo motions are modeled by motions 
of density peaks of a suitably smoothed version of the ini- 
tial Gaussian density field evolved according to linear theory 
(Peebles 1980). iV-body simulations have been used to test 
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those predictions (e.g. Bahcall, Gramann, Cen 1994; Croft 
& Efstathiou 1994; Suhhonenko & Gramann 1999, Colberg 
et al. 2000; Sheth & Diaferio 2001; Yoshida, Sheth & Diafe- 
rio 2001). Colberg et al. (2000) found that, at z > 4, the 
mean peculiar velocities of materials that end up in mas- 
sive halos agree well with the predictions of linear theory 
for the corresponding peaks, but that linear theory system- 
atically underestimates their subsequent growth; the rms 
peculiar velocities at z = evaluated from TV-body simu- 
lations are larger than the linear theory prediction by up 
to 40 percent. This discrepancy may be ascribed to the fol- 
lowing facts: that evolution of peculiar velocities of halos in 
high density regions are accelerated relative to the average 
(Tormen & Bertshinger 1996; Cole 1997; Kepner, Summers 
& Strauss 1997) and that halos tend to be formed in such 
higher-density regions because of the biased structure forma- 
tion mechanism. Sheth & Diaferio (2001) have shown using 
TV-body simulations that the evolution of the halo peculiar 
velocity indeed depends on the local matter density, rather 
than on the global density. They presented a simple model 
for the shape of the distribution function of the peculiar ve- 
locities that takes into account their dependence on the local 
density in an empirical manner. 

The purpose of this paper is to present a simple and 
accurate model for the rms of peculiar velocities of halos. 
Hereafter, we denote the rms peculiar velocity of halos with 
mass M as crhaio(M). 

The major improvements that we achieve here are to 
take into account the following two points that play an im- 
portant role in understanding statistical properties of halo 
peculiar velocities, (i) The evolution of the halo peculiar ve- 
locity depends on the local matter density, rather than on 
the global density. We develop a model of the local density- 
dependent growth of the halo peculiar velocity in an em- 
pirical manner following Sheth & Diaferio (2001). (ii) Dark 
matter halos are biased tracers of the dark matter distri- 
bution, and thus halos tend to be located at high-density 
regions. In order to implement this, we use the halo biasing 
model developed by Mo & White (1986) and Sheth & Tor- 
men (1999; 2002) combined with an empirical log-normal 
probability distribution function (PDF) of the mass density 
distribution (Kayo, Taruya & Suto; 2001). 

The outline of this paper is as follows. In Section 2, we 
summarize linear theory predictions for the peculiar veloci- 
ties of peaks. In Section 3, we describe TV-body simulations 
that are used to calibrate and to test the model. We con- 
struct the model in Section 4. The model predictions are 
compared with TV-body simulations in section 5. A summary 
and discussion are presented in Section 6. 
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The linear growth factor D(a) is given by 
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Note that the frequently used approximation f(a) = 
d In D/d In a ~ Q, 0,6 underestimates the growth rate by ~ 5 
percent in the flat fio = 0.3 cosmology that we consider here. 
Thus, we will not adopt the approximation but we will use 
the numerical integration throughout the present paper. 

We assume that the primordial (linear) density fluctua- 
tions obey the Gaussian statistics. Then the peculiar veloc- 
ity field is isotropic and Gaussian, and its three-dimensional 
dispersion smoothed over a smoothing scale of ii is given by 



a v (R, a) = df(a)a-i(R, a). 



(4) 



Following Bardeen et al. (1986), Oj is denned for any integer 
j as 



1 

2^2 



dk k 2+23 P(k)W z {kR), 



(5) 



where W(x) is the Fourier transform of the smoothing 
window, and P(k) is the linear matter power spectrum. 
Throughout this paper, we use the real space top-hat win- 
dow function which corresponds to W{x) = (3/a; 3 )[sin(x) — 
a;cos(a;)]. We adopt the fitting function of a cold dark mat- 
ter (CDM) power spectrum given by Bardeen et al. (1986), 
who show that the rms peculiar velocity at the peaks of the 
smoothed density field differs systematically from <r„ by 



a p (R) = a v {R)^l-^/ala\. 



(6) 



Note that this expression does not depend on the height of 
peaks. 

The conventional assumptions in modeling the peculiar 
velocities of halos are as follows, (i) The velocities of halos 
are identical to those of peaks (equation The smoothing 
scale is chosen such that the filtered mass is equal to that 
of the halo considered, i.e., M = 4,-KpR? /i in the case of 
the real space top-hat window function, (ii) The halo veloci- 
ties evolve according to linear theory (equation^. However, 
the prediction based on these assumptions seems to under- 
estimate the rms peculiar velocities of halos (Colberg et al. 
2000). 



2 LINEAR PREDICTIONS FOR THE 
PECULIAR VELOCITIES OF PEAKS 

In this section, we summarize basic equations of the linear 
theory predictions for peculiar velocities of density peaks. 

According to linear theory of gravitational instability in 
a dust universe (Peebles 1980) , the peculiar velocity field of 
mass grows as 

v oc aV = au — 

a In a 



3 TV-BODY SIMULATION 

In constructing an improved model of the peculiar veloc- 
ity of halos, we first use TV-body simulations of a ACDM 
cosmology (Qq = 0.3, Ao = 0.7, h — 0.7, and as — 1.0). 
The details of the TV-body simulations are described in Jing 
& Suto (1998); see also Kayo et al. (2001). Briefly, the 
simulations employ 256 3 CDM particles in a cubic box of 
300/i _1 Mpc on a side. The simulations are performed us- 
ing a P 3 M code with the gravitational softening length of 
e ~ 120 /i kpc. The initial matter power spectrum is com- 
puted using the CDM fitting function given by Bardeen et 
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Figure 1. Halo velocity dispersions as a function of halo mass 
at different redshifts as indicated in each panel. Filled circles are 
measurements from JV-body simulations. The long-dashed lines 
represent the phenomcnological model, equation 10 adopting the 
conditional mass function via equation J5J, while solid lines are 
for the model with the approximation equation ill 01 . The dashed 
lines show the linear theory prediction for peaks (equation [U . 

al. (1986). We have three realizations which differ only in the 
phases of the initial density fluctuations. The particle mass 
(rripart = 1.34x h~ x Mq) of the simulations is sufficiently 
small to guarantee that there are practically no discreteness 
effects on dark matter clustering on scales down to the soft- 
ening length in the redshift range of interest for our purposes 
(Hamana, Yoshida & Suto 2002). 

We identify dark matter halos using the standard 
friends-of-friends algorithm with a linking parameter of 6 = 
0.164 (in units of the mean particle separation). We set the 
minimum mass of the halos as 1.34x 10 12 /i -1 Mq, which cor- 
responds to the mass of 10 simulation particles. We define 
the peculiar velocity of each halo to be the center-of-mass 
peculiar velocity of member particles. 

The smoothed mass density fields are computed on 64 3 
regular grids in the simulation box using the top-hat window 
function. For each halo, the smoothed mass density of the 
nearest grid is assigned as its local background density. 

Fig. shows the rms peculiar velocities of halos in the 
A-body simulation (filled circles) together with the linear 
theory prediction of peaks, (equation |S] dashed curves) as 
well as our predictions described in the next section (solid 
lines and long-dashed lines). Note that, in order to account 
for the finite size of the simulation box (Lbox = 300/i -1 Mpc), 
we set a lower limit of the integration in equation JSJ by 
k m in = 27r/Lbox- This large-scale cut-off reduces a p by a 
fractional amount of about 8 percent, and is applied to all 
model predictions in this paper. The plot confirms the previ- 
ous finding that the linear theory predictions underestimate 
the rms peculiar velocities by about 10-40 percent for cases 
< z < 2.2 in a ACDM cosmology. In addition, it is seen 
from this Figure that the linear theory predictions have a 



slightly steeper slope than that of A-body simulations, es- 
pecially for the case of z = 0. This discrepancy may be ex- 
plained by the following: (i) massive halos tend to be formed 
in denser regions as a result of the biasing mechanism, and 
(ii) halos in denser regions move faster than those in less 
dense regions, as will be proved in the next section (Fig. 
and^J. As a consequence, on average, deviations from linear 
theory appear larger for massive halos, and thus the slope 
of A-body data becomes flatter. 



4 MODEL 

In this section, we describe our model of the rms peculiar ve- 
locity of halos. We basically follow the framework developed 
by Sheth & Diaferio (2001); see also Diaferio & Geller (1996) 
and Sheth (1996). They originally developed the model for 
the distribution function of peculiar velocities of dark matter 
particles and halos, but their framework can be straightfor- 
wardly applied to the model for the rms of halo peculiar 
velocities. 

Sheth & Diaferio (2001) pointed out that the evolution 
of the halo peculiar velocity depends on the local matter den- 
sity (defined with an appropriate smoothing length which we 
discuss later) . We thus first construct a model for the depen- 
dence of the halo peculiar velocity dispersion, Oh alo (M, S), on 
the local background density 8. Following Sheth & Diaferio 
(2001), we consider a parametric model 

oLo(M, 5) = [1 + S(R local )] 2 » (R '°^a 2 p (M), (7) 

where i?i oca i is the smoothing scale with which the local 
background density S is defined. The key question is how 
to define the appropriate smoothing length-scale. Evidently, 
-Riocai should enclose the gravitational coherence scale which 
is responsible for the local deviation of the peculiar motion 
of a halo from its global value (for instance, given by linear 
theory). Thus, the size of halo (e.g., the virial radius) sets 
a lower limit for 7?i oca i. On the other hand, choosing values 
for -Ri oca i that are too large would result in an excessively 
smoothed density field, and local density fluctuations, by 
which the deviation from the global value is induced, are 
smoothed out. We may, thus, expect that i?i OC ai is in the 
range between the linear scale (such that ao(R) < 1) and 
the non-linear scale (ao(R) > 1)- Concerning the power- 
law index /i, the linear theory relation, / cx f2 0,6 , suggests 
[i ~ 0.6 (Sheth & Diaferio 2001). Note, however, that fj, 
depends on 7?i OC ai on which the local density is defined (a 
larger /i for a larger smoothing length, because the range of 

5 becomes narrower), thus \i does not need to be very close 
H ~ 0.6, but it is expected that with an appropriate choice 
of the smoothing length-scale, /i would be close to 0.6. 

To proceed further, we adopt an ansatz that i?i oca i is 
given by a relation cro(-Riocai) = o"i oca i, and we determine 
the model parameter <ri oca i empirically using A-body sim- 
ulations. As expected from the discussion above, a suitable 
choice of ai oca i would be within a range 0.1 < (Ji OC ai < 1. We 
find that choosing /i = 0.6 (constant) with i?i oca i given by 
0-iocai = 0.5(1 + z)~ - 5 (through cro(-Riocai) = criocai) provides 
reasonable fits to the results of the A-body simulations as 
shown in Fig.|2] For this choice, Ri oc &x — 18 and 16/i _1 Mpc 
for z=0 and 0.8, respectively. The value of /x = 0.6 is con- 
sistent with one suggested in linear theory (/x ~ 0.6). Thus, 
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Figure 2. Halo peculiar velocity dispersions as a function of 
the local background density for limited mass ranges. Halo mass 
ranges (in units of h~ Mq) are denoted in each plot. Dotted lines 
represent the model equation J7J with fi = 0.6 and filled circles 
show measurements from our TV-body simulations. The smoothed 
local density fields are computed from the TV-body particle dis- 
tributions using the top-hat window function with a smoothing 
scale of 20ft, -1 Mpc. 
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Figure 3. Same as Fig. l2*lbut for halo progenitors at the initial 
redshift of the TV-body simulation (z = 18). We mark the member 
particles of a halo identified at z = 0, and then compute the 
position and velocity of the halo progenitor from the center-of- 
mass position and velocity of the marked particles at z = 18. 
The top-hat smoothing scales for computing the local density are 
taken to be 20/i _1 Mpc. 



while the choice of the smoothing length may still be ar- 
bitrary, we decide to use the above sets of parameters for 
definiteness. 

Fig.[3]shows the peculiar velocity dispersion of halo pro- 
genitors as a function of the local density at the initial red- 
shift of the TV-body simulation z = 18. We mark the mem- 
ber particles of a halo identified at z = 0, and then compute 
the position and velocity of the halo progenitor from the 
center-of-mass position and velocity of the marked particles 
at z = 18. Fig. [3] clearly indicates that there is no clear 
correlation between the peculiar velocity and the local den- 
sity at the initial epoch. Together with Fig.|3 this strongly 
supports the idea that most of the halo peculiar velocity 
is acquired through the subsequent evolution of the nearby 
density field. A closer look at Fig. |31 reveals that the initial 
peculiar velocity dispersions decline very slightly at both the 
lowest and highest bins. An interpretation of this is that, at 
the initial epoch, the halo progenitors locate close to posi- 
tive (negative) density peaks tend to have relatively higher 
(lower) local densities. At the same time, a potential gradi- 
ent at positions close to the density peaks tends to be small, 
and thus the progenitors at such regions move slower than 
others. As a result, the highest and lowest density bins have 
a slightly lower peculiar velocity than other bins. 

As shown by Sheth & Diaferio (2001), if the range 
of halo masses and local background density is sufficiently 
small, the PDF of halo peculiar velocities is well approxi- 
mated by a Maxwellian distribution, i.e., the PDF of one- 
dimensional velocity is a Gaussian. We denote by p(M\S) 
the probability of finding a halo with mass M at a region 
where the background density is S. Then the peculiar veloc- 
ity dispersion of halos with mass M is given by summing up 
the dispersion cr^ alo (M, S) weighted over the probability of 
finding halos in regions of 8: 



Ohalo(M) 



J d8 p{M\8)oi^{M,5) 
J d8 p(M\8) 



(8) 



Note that this expression is valid only if the background 
density is defined appropriately (i.e., with the smoothing 
scale -Riocai defined above). The expression for p(M\8) is 
obtained as follows. 

It is well established that dark matter halos are biased 
tracers of the underlying dark matter distribution. In this 
picture, it is useful to write the conditional probability as 



(9) 



where n(M) and n(M\8) denote the halo mass function and 
the conditional mass function in regions with the local back- 
ground density 8, respectively, and Pdm{8) is the PDF of the 
dark matter density field. To evaluate the first term of the 
right hand side of equation ©, i.e. the excess number den- 
sity of halos, we adopt two approaches. The first is to directly 
use the mass functions derived from the extended model of 
Press-Schechter (1974); see, for example, Bond et al. (1991). 
Lacey & Cole (1991) and Mo & White (1996). In practice, 
we adopt the conditional mass function given by Sheth & 
Tormen (2002) that even takes into account the effect of the 
ellipsoidal collapse. The other approach is to use a simple 
approximation, n(M\S) ~ [1 + b(M)S]n(M). Note that this 
approximation is valid only if the smoothing scale defining 
the local density is much larger than the halo size, or equiv- 
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Figure 4. Conditional mass functions normalized by the uncon- 
ditional mass function for ranges of the local background density, 
from lower to upper, — 1 < 8 < 0, < 8 < 1, 1 < d < 2 and 
2 < 8 < 3. Dashed curves show the conditional mass function by 
Sheth & Tormen (2002), dotted curves show the approximation 
n(M\8) ~ [1 + b(M)8]n(M) with the unconditional mass function 
and halo bias model by Sheth & Tormen (1999), and the solid his- 
tograms are for ./V-body data. The top two cases are offset upward 
by a factor of 10 and 100. The smoothing scales for computing 
the local density are -f?i oca i = 20 and 10h~ x Mpc for z = (left 
panel), and 0.8 (right panel), respectively. 



alently if the total mass enclosed within the smoothing scale 
is much larger than the halo mass. Adopting this, equation 
10 reduces to 

p{M\5) = [1 + b(M)5] Pdm {5). (10) 

Mo & White (1996) developed a Lagrangian halo bias model 
based on the extended Press-Schechter model. Jing (1998) 
and Sheth & Tormen (1999) discussed the correction for the 
mass dependence of the halo bias and proposed modified 
fitting functions that give better fits to numerical simula- 
tion results. In what follows, we use the linear bias model 
of Sheth & Tormen (1999). Although the former approach 
may lead to a better model than the latter, we examine the 
latter approximation to see if the simple approach could 
provide an acceptable model. Fig. [I] compares the above 
two theoretical predictions (dashed and dotted lines) for 
n{M \S)/n(M) with the A r -body data (solid histograms) for 
ranges of the local background density. The simple approx- 
imation n(M\S) ~ [1 + b{M)S\n{M) turns out to provide 
a good fit to the iV-body data for smaller mass ranges 
Mh a i < lO 13 ~ 14 ft -1 M0, but underestimates the conditional 
probability for larger mass ranges. The conditional mass 
function by Sheth & Tormen (2002) gives good fits for high 
local densities (8 > 1, upper two cases in the plots), while 
fits become degraded for low densities (5 < 1, lower two 
cases in the plots). Importantly, it reproduces the cut-off 
in the conditional mass function at a large mass nicely in 



marked contrast to the approximation which does not have 
the cut-off. 

The dark matter PDF is known to be well fitted to the 
log-normal distribution (e.g., Coles & Jones 1991; Kofman 
et al. 1994; Kayo, Taruya & Suto 2001; see also Taruya, 
Hamana & Kayo 2003 for a possible explanation of the origin 
of the log-normal PDF) 

/« 1 / [ln(l + (5)+a?/2] 2 ) 1 

Pln{5) = j2^y7t e * p \ 2o{ jiTs 

where <ti is related to the variance of the nonlinear density 
field by 

cr*(R) = ln[l + a 2 nl (R)}, (12) 
with 

°ni(R) =2^f dk k 2 P n i(k)W 2 (kR). (13) 

In the above expression, P n i(k) denotes the nonlinear mat- 
ter power spectrum for which we compute using the fitting 
function of Peacock & Dodds (1996). 



5 RESULTS 

Fig. U compares our model predictions with TV-body simula- 
tions (filled circles) ; long-dashed curves correspond to equa- 
tion © with the conditional mass function by Sheth & Tor- 
men (2002), while solid curves correspond to the approxima- 
tion equation I1U1 . It is clear that our models (especially the 
former one) substantially improve the accuracy of the pre- 
diction compared with that of linear theory (dotted lines) 
except for a very high mass range of M > 10 14 h^ 1 Mq. How- 
ever, the predictions of the model adopting the approxima- 
tion equation ffllUB rise very rapidly at a high mass range of 
M > 10 14 Mq and significantly deviate from the simula- 
tion results. This is explained as follows: As shown in Fig. [I] 
the conditional mass function computed from A r -body sim- 
ulation declines at a high mass range. However, the approx- 
imation equation 1 1 (H does not reproduce this feature but 
keeps on rising toward a large halo mass. Consequently, the 
approximated model over-estimates the contribution from 
moderately over density regions. Here, it must be noted that 
with our choice of the smoothing scale (i?iocal = 18, 15 and 
13/i _1 Mpc for z — 1, 2 and 3, respectively), the most of 
the matter in the universe is in regions where the local den- 
sity contrast being — 1 < S < 1. Therefore, a discrepancy 
in the conditional mass function at a high density range of 
5 > 1 between the model and the simulation does not have 
a large influence on the model prediction. Adopting the con- 
ditional mass function proposed by Sheth & Tormen (2002) , 
which reasonably reproduces the decline as shown in Fig. [I] 
improves the accuracy of the prediction very well. In par- 
ticular, this model predicts a flatter slope than that of the 
linear theory prediction, and it is very close to the simula- 
tion results. We, however, note that the model prediction 
adopting Sheth & Tormen conditional mass function rises 
at very large mass scales (M > 10 15 ) and deviates from the 
simulation. The reason of this is that the conditional mass 
function becomes less accurate at such mass range. Because 
at that mass range, the number density of halos is too small 
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Figure 5. PDFs of one-dimensional halo velocity v\ at z = (left 
panels) and z = 0.8 (right panels) and for limited mass ranges 
(denoted in each plot). Histograms show the measurement from 
the TV-body simulations. Solid curves show our model prediction, 
while dotted curves represent the Gaussian distribution with a 
computed from the TV-body simulation. 
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Figure 6. Pairwise velocity dispersion of halos as a function of 
the separation for limited mass ranges (denoted in each plot). 
Filled circles indicate the TV-body results. Dotted lines represent 
the linear theory prediction for peaks neglecting velocity correla- 
tion, {v^) 1 ^ 2 = ^/2/3(7p , while solid lines show our model pre- 
diction («i 2 ) 1/2 = \/273o- halo . 



to have a meaningful statistic, this disagreement is not seri- 
ous. 

Although there is still room for improvement in the 
model, our model adopting the conditional mass function 
of Sheth & Tormen (2002), provides a good fit to the N- 
body data over a wide range of the halo mass. To improve 
the model prediction, it may be necessary to adopt a fur- 
ther developed conditional mass function. Developing an ac- 
curate conditional mass function is beyond the scope of the 
present paper, and hence we do not attempt this here. In the 
following discussions, we adopt the model adopting Sheth & 
Tormen conditional mass function. 

As Sheth & Diaferio (2001) originally showed, the model 
presented in this paper is applied to model the PDF of halo 
peculiar velocities 



P(v\M) = 



J d5 p(M\5)P(v\M,5) 
JdS p(M\S) ' 



(14) 



where P(v\M) is the PDF of peculiar velocities of halos 
with mass M, and P(v\M,S) is the PDF for halos at re- 
gions with the local background density 8. The PDF of 
one-dimensional peculiar velocities (which we denote by Vi) , 
P(vi\M, 8), is well approximated by the Gaussian distribu- 
tion (Sheth & Diaferio 2001). We compute a of the Gaus- 
sian PDF using equation (J7) with a reasonable isotropic as- 
sumption (v 2 ) = (v 2 )/3. The results are plotted in Fig. 
(solid lines), and are compared with TV-body simulations 
(histograms) and Gaussian PDF (dotted lines). As Fig. |S] 
clearly shows, our model agrees with the TV-body simula- 
tion very well. Also it is seen from the Figure that the PDF 
of one- dimensional peculiar velocity consists of two parts; a 
Gaussian core and an exponential tail. The exponential tail 
arises from the sum of Gaussian distributions with differ- 



ent dispersions (Sheth & Diaferio 2001; Kuwabara, Taruya 
& Suto 2002). Our model nicely reproduces the deviation 
from the Gaussian distribution at large peculiar velocities 
and agrees with the simulation very well. 

The model presented in this paper can be also applied 
to the pairwise velocity dispersion of halos at large dis- 
tances. The model of the pairwise velocity dispersion plays 
an important role in computing the redshift distortion of 
the two-point correlation function of halos. Because any ve- 
locity correlation between halos with large separations of 
r > 20/i _1 Mpc may be safely ignored, the halo pairwise ve- 
locity dispersion is given by {v 2 2 ) = 2/3 (v 2 ), where viz de- 
notes the pairwise velocity of a pair of halos. Often a 2 (equa- 
tion^ is used to compute {v 2 ) (e.g., Hamana et al. 2001); 
see Sheth et al. 2001a,b for an alternative approach. How- 
ever, as shown in the previous sections, the linear theory 
prediction underestimates the halo velocity and our model 
provides better fits to TV-body data. We may, therefore, ex- 
pect that our model improves predictions of the pairwise ve- 
locity dispersion. Fig. compares the model prediction for 
the pairwise velocity dispersion (solid lines) with the result 
of our TV-body simulation (filled circles). We also plot the 
linear theory prediction for peaks (dotted lines) . As Fig. 
clearly shows, our model indeed improves the accuracy of 
prediction at r > 20/i _1 Mpc. 



6 SUMMARY AND DISCUSSIONS 

We present a model of peculiar velocities of dark matter 
halos which describes various statistical properties well and 
significantly improves the accuracy of the simple linear the- 
ory predictions, including the halo rms velocity dispersion, 
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PDF, and the pairwise velocity dispersion at large separa- 
tions. The model predictions, however, becomes less accu- 
rate for very large halo masses (M > 10 15 ft -1 Mq). This 
is mainly due to the limitation of the adopted conditional 
halo mass function, and may be reconciled by using an im- 
proved model. Note that this shortcoming hardly causes any 
practical problem, because the number density of such very 
massive halos is too small to have a meaningful statistic. 

The model presented in this paper is characterized by a 
single free parameter, /x(-Riocai), that we determine in an em- 
pirical manner. To do this, we need the appropriate smooth- 
ing scale for the current purpose. We make an anzats that 
J^iocai is in the range between the linear and the non-linear 
scale and it should be given by the relation (Xo(-Rlocal) = 
ciocai for some value in a range 0.1 < ai oca i < 1. This turns 
out to be reasonably successful, and we have found that a 
choice of /i = 0.6 with o"i oca i = 0.5(1 + z)" ' 5 works well for 
all redshifts we consider in this paper (z < 2.2). This choice 
gives the smoothing scales of -Ri oca i = 18, 15, 11, 8.5 and 
5.6/i _1 Mpc for z = 0, 1, 3, 5 and 10, respectively. However, 
we calibrated the model with ./V-body data at < z < 2.2 
only, thus the validity of this procedure is questionable for 
very higher redshifts (although in practice velocities of ha- 
los at very high redshifts are not observable currently). A 
more sophisticated approach to the choice of the smoothing 
scale is needed for a further improvement of the model. This 
may be closely connected to an understanding of the role of 
the local environment in the structure formation, and is a 
subject of work in progress. 
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